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We have calculated the contributions proportional to /3q a™ +1 to the longitudinal fragmentation function in 
e + e~ annihilation to all orders of perturbation theory. We use this result to estimate higher-order perturbative 
corrections and nonperturbative corrections to the longitudinal cross section ox and discuss the prospects of 
determining a a from ox. The structure of infrared renormalons in the perturbative expansion suggests that 
the longitudinal cross section for hadron production with fixed momentum fraction x receives nonperturbative 
contributions of order l/(x 2 Q 2 ), whereas the total cross section has a larger, 1/Q correction. This correction 
arises from very large longitudinal distances and is related to the behaviour of the Borel integral for the cross 
section with fixed x at large values of the Borel parameter. 



1. Introduction 

The ALEPH § and OPAL § collaborations 
have measured the dependence of single-particle 
inclusive cross sections in e + e~ annihilation on 
the scattering angle 8 between the observed 
hadron h and the incoming electron beam. The 
angular dependence discriminates between contri- 
butions from transversely and longitudinally po- 
larized virtual bosons, and from Z°-photon inter- 
ference S 
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In the following, dropping the superscript '/i' im- 
plies summation over all hadrons h. 

In this paper we concentrate on the longitudi- 
nal cross section. It is given as a convolution of 
a parton fragmentation function (p = q, q, g) 
with a partonic cross section da v L j dx 
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The perturbative expansion of the longitudinal 
parton cross section starts at order a s . 

Summed over all hadrons, the fragmentation 
functions satisfy the energy conservation sum rule 
J^dx xDp(x,Q 2 ) — 1. Consequently, the in- 
tegrated longitudinal cross section is an infrared 
(IR) safe quantity which is calculable in pertur- 
bation theory 
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Here <7q is the Born total e + e" annihilation cross 
section, a s = a s (Q) and Nf is the number of ac- 
tive fermion flavours. The next-to- leading order 
contribution has been obtained in Q). OPAL Q 
has measured <jl the Zq peak: 

<J L /<Jtot(M 2 z ) = 0.057 ± 0.005 . (4) 

One of the main motivations for the present 
study is to investigate whether measurements of 
the total longitudinal cross sections can yield 
a precise determination of the strong coupling. 
This requires that we control higher-order per- 
turbative corrections and nonperturbative effects, 
both of which are expected to be much larger 
for ctl than for the total cross section a to t- We 
address both types of corrections in this report, 
by studying the structure of IR renormalons, a 
certain class of higher-order perturbative correc- 
tions, for the longitudinal cross section. 

The study of nonperturbative effects in frag- 
mentation functions is an interesting topic in its 
own right . The light-cone expansions for frag- 
mentation functions and for structure functions in 
DIS are similar and suggest that nonpertur- 
bative effects in both cases are of order 1/Q 2 and 
can be described in terms of multi-parton distri- 
butions. In contrast to DIS, however, the relevant 
operator structures for fragmentation are essen- 
tially nonlocal and cannot be expanded at small 
distances. Hence the usual operator product ex- 
pansion does not apply and the status of the light- 
cone expansion is less established. Hadroniza- 
tion models generically introduce nonperturba- 
tive corrections of order 1/Q, while current data 
on scaling violations in fragmentation do not dis- 
tinguish between 1/Q or 1/Q 2 behaviour. A non- 
perturbative correction of order 1/Q to the to- 
tal longitudinal cross section was suggested in 
§ as a consequence of phase-space reduction in 
the one-loop diagram calculated with a massive 
gluon. We address these apparently conflicting 
statements below. 

While this work was in writing, Dasgupta and 
Webber ||] have addressed a similar set of ques- 
tions with closely related methods. 




Figure 1. The 'primary' quark contribution to 
ul ■ Sum over all possible insertions of the bubble 
chain is understood. 



2. General formalism 

There is suggestive evidence from exact low- 
order results that /3o is a large parameter and that 
keeping corrections of order ([3oa s ) n in higher 
orders resums important contributions. More- 
over, the infrared renormalons encoded in the 
corresponding series can elucidate the power- 
behaviour of nonperturbative corrections and, 
perhaps, even their x-dependence, as in the case 
oida^/dx. The {(3oa s ) n corrections can be traced 
by inserting a chain of fermion loops into the 
gluon propagator, and by restoring the full QCD 
/3-function coefficient /3 = -1/(4tt)[11 - 2/3 N f ] 
from the dependence on Nf. For ul we obtain 
two contributions, according to whether the reg- 
istered parton comes from the primary vertex or 
a fermion loop. The corresponding diagrams are 
shown in Fig. 1 and Fig. 2, for the contributions 
of the 'primary' and 'secondary' quarks, respec- 
tively. Note that the secondary quark contribu- 
tion reduces to the gluon contribution at lowest 
order in a s . 

The evaluation of the two classes of diagrams, 
for an arbitrary number of internal fermion loops, 
and for their sum, is relatively straightforward by 
means of the dispersion technique developed in 
fl0|-|l2|| , in terms of the distribution function over 
the invariant mass k 2 of the bubble chain: 

j~[P,s] 

^( X ^ = k 2 /Q 2 )^ 
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Figure 2. The 'secondary' quark contribution to 
<jl- Sum over all possible insertions of the bubble 
chain is understood. 
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The notation for momenta and Lorenz indices 
corresponds to Fig. 2. We denote by [p] ([s]) 
the contribution of the 'primary' ('secondary') 
quark, while JdLips[. ..] are Lorentz-invariant 
phase space integrals with the momentum con- 
servation 5-function included, and M^-M*,^ is 
the matrix element for the primary qqg amplitude 
squared. 

Note that in the case of the 'primary' quark 
contribution the phase space integral over k\ , ki 
is proportional to k 2 , so that the result takes 
the form of the one-loop diagram calculated with 
a gluon of mass k 2 . This equivalence does not 
hold for the registered 'secondary' quark because 
of the nontrivial longitudinal projector. Because 
of this inequivalence, the restoration of (3q from 
fermion loops is not unambiguous and the relation 
of fermion loop chains with running coupling ef- 
fects is partially lost. In practice, we have found 
that the numerical differences arc small, so that a 
detailed discussion is deferred to H . The analytic 
expressions for the primary and secondary quark 



contribution are rather lengthy and will also be 
given there. 

Given the invariant mass distributions (in £ = 
fc 2 /Q 2 ), finite order results are obtained in terms 
of the logarithmic integrals jlO 12 
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The sum of the series, defined by a principal value 
prescription for the Borel integral, equals 
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where £l < is the position of the Landau pole in 
the strong coupling and the function $(£;) is spec- 
ified in Eq. (2.25) of |Q. Infrared renormalons 
correspond to nonanalytic terms in the expansion 
of da l I 'dx{x , £) at small £ 
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and are interpreted as indications of nonpertur- 
bative power corrections of the form 
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Their size can be estimated by the corresponding 
ambiguity in the summation of the perturbative 
series, which is of order of the imaginary part (di- 
vided by 7r) of the sum in (]?]). Note that identify- 
ing the x-dependence of the power corrections in 
(||) with the x-dependence of the IR renormalon 
ambiguity or, equivalently, the coefficients of non- 
analytic terms in (^) is an assumption which can 
not be justified from first principles. Since IR 
renormalons in short- distance quantities are re- 
lated to ultraviolet ambiguities in higher-twist 
matrix elements, we refer to this assumption as 
the 'ultraviolet dominance' of higher-twist correc- 
tions. 
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3. Perturbative series for 07, 

In this section we consider perturbative correc- 
tions to <Tl, written as 



(10) 



where 00 is the Born total e + e~ cross section. 
As mentioned earlier, we approximate the exact 
higher-order coefficient by its value in the 'large- 
/?o' limit, where /?o is restored from the term with 
the largest power of Nf at each order. This ap- 
proximation, called 'naive nonabelianization' in 
HQ] i reduces to the familiar BLM prescription for 
n=l. To see how it works, we rewrite the exact 
a 2 , correction in (^|) as 

dx =6.17 - 0.7573/(-A>). (11) 

With — /3o = 0.61 for Nf — 5, neglecting the sec- 
ond term gives an accuracy of about 25%. We 
have calculated the coefficients d n in higher or- 
ders, in the MS scheme. The 'primary' and 'sec- 
ondary' quark contributions, d$ and dn , respec- 
tively, add to dn as d n = dm /3 + 2dn /3. A few 
lower order results up to order a* are 
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The sum of these contributions to all orders is 
conveniently written in terms of 'enhancement 
factors' relative to the leading order contribution 
111 defined by 



M^(a s ) = l + Y,(-Poa s ) n dt 
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For various values of a s (Mz) we get at Q = Mz 
a s = 0.110 : 

M [s] = 1.92 ±0.05 . 



MM = 1.59 
a s = 0.120 : 

MM = 1.68 
a s = 0.130 : 

MM = 1.79 



M [s] = 2.08 ± 0.08 
M [s] = 2.23 ±0.12 
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Figure 3. Longitudinal fraction in the total 
e + e _ cross section: (c) leading order, (b) next-to- 
leading order and (a) resummation of all orders in 
/?qQ!™ +1 corrected for the exact 0(a 2 s ) coefficient. 



The given numbers correspond to a principal 
value definition of the Borel integral and the un- 
certainties roughly coincide with the size of the 
minimal term in the series^. Let us add the fol- 
lowing comments: 

(i) The perturbative coefficients in (|l2[), ( |l3| ) 
grow rapidly, especially for the secondary quark 
contribution. This growth is related to an IR 
renormalon, that indicates a l/Q 2 correction to 
primary quark fragmentation and a l/Q cor- 
rection to secondary quark fragmentation, see 
Sect. 4. 

(ii) Even though the l/Q power behaviour in- 
dicates much larger nonperturbative corrections 
to <tl as compared to er tot , the moderate size of 
the minimal term of the perturbative series sug- 
gests that these corrections are still not large at 
Q = Mz- The relatively large hadronization cor- 
rection for <jl within the JETSET model applied 
in could thus correspond to higher-order per- 
turbative rather than nonperturbative effects. 

(iii) This suggestion is also supported by Fig. 3, 
where for a s (Mz) = 0.118 we have plotted the 
energy dependence of the total longitudinal cross 
section. Taking into account higher-order pertur- 
bative corrections [curve (a)] steepens the energy 
dependence, such that it is not far from the JET- 
SET prediction, where the steep energy depen- 



(16) 



4 The corresponding uncertainty for M ' p l is small in com- 
parison with the one for MW and is omitted. 
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dence is due to the hadronization correction. It is 
worth noting that the parton shower Monte Carlo 
alone does not yield this energy dependence. Ex- 
perience with similar calculations suggests that 
the approximation of resumming only (f3oa s ) n 
contributions overestimates radiative corrections, 
so that we expect a more realistic estimate in be- 
tween the curves (a) and (b). An exact O(a^) cal- 
culation would reduce the theoretical error con- 
siderably. 

(iv) In universality of the l/Q-power cor- 
rection was assumed and a corresponding unique 
phenomenological parameter fitted from the dif- 
ference between the measured average thrust (1 — 
T) and the theoretical second order prediction. 
When added to the second order result (||), one 
obtains a prediction for ctl consistent with data. 
There is no conflict between the procedure of [ [L3[ 
and the one presented here, if the phenomenolog- 
ical 1/Q correction effectively parameterizes the 
higher-order perturbative contributions added in 
our approach. If universality of power corrections 
holds, these perturbative corrections would also 
be universal, at least asymptotically in large or- 
ders. However, from the point of view presented 
here, the universality assumption is not required, 
since higher-order corrections are in principle cal- 
culable for each observable. 

4. Power corrections 

Returning to we quote the expansions for 
the invariant mass distributions in (||): 
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Interpreting £ as (A/Q) 2 where A is the QCD 
scale, these expressions are valid for x > A/Q. 



We note that for such x, all power corrections 
are at most of order 1/Q 2 , in agreement with 
the result from the light-cone expansion of frag- 
mentation processes in [|J. We also see that 
the power expansion runs in A 2 /(Q 2 x) for the 
primary quark contribution and A 2 / \Q 2 x 2 ) for 
the secondary quark (gluon) contribution. The 
strong divergence of the second contribution for 
small x makes it possible for the moments of 
the x-distribution to have parametrically larger 
power corrections. Indeed, we find for the two 
contributions to the total longitudinal cross sec- 
tion 
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with a 1/Q correction for the secondary quark 
contribution. Assuming ultraviolet dominance 
of higher-twist corrections, the x-distributions 
given in (p"7j), JT8] ) can be used to model the x- 
dependence of power corrections by convoluting 
the partonic power correction with the leading 
twist fragmentation function MM. Note that the 
expressions for the secondary quark contribution 
differs from the gluon contribution to <tl in 
because the series of higher-order fermion loop di- 
agrams does not reduce to the massive gluon cal- 
culation performed in j^] . The ensuing additional 
model dependence in the estimate of higher-twist 
corrections will be discussed in |^| . Both the cal- 
culation here and the calculation with a massive 
gluon coincide in the essential aspects — power 
corrections of order A 2 /(Q 2 x 2 ) for finite x and 
1/Q for the integrated longitudinal cross section. 

In the remainder of this section, we discuss 
the origin of the 1/Q correction in more detail. 
Adopting for this purpose the massive gluon ap- 
proximation, one finds that the Borel transform 
for the gluon fragmentation contribution factor- 
izes as 
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when some terms that can not give rise to a 1/Q 
correction are omitted. Here u is the Borel pa- 



G 



rameter and analyticity of F(u) for |u| < 1 cor- 
responds to the statement that only 1 /Q 2 correc- 
tions arise at finite x. Now, the ^-dependence can 
be absorbed completely into a change of scale in 
the coupling and the Borel integral is given by 



takes the form (a = 1 — a) 



do 
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du exp 
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The Borel integral (leaving renormalon poles at fi- 
nite u aside) does not exist for x < A/Q, because 
it diverges at infinity. This is a manifestation 
of the fact that for such small x the power ex- 
pansion breaks down and that power corrections 
to integrated distributions depend sensitively on 
how the small- x region is weighted. Indeed, be- 
cause of the factorization of the x-dependence, 
the ^-integration is trivial and we get 
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For the total longitudinal cross section, 7 = 0, 
and the newly generated pole at u = 1/2 cor- 
responds to the 1/Q correction discussed before. 
Note that effects due to color coherence and an- 
gular ordering are expected to change the small- 
x asymptotic behaviour, which could potentially 
shift the pole to a different value. Clarifying 
the impact of resummation of small- x logarithms 
requires similar efforts to those that have been 
undertaken to understand the effect of Sudakov 
resummation on power corrections in Drell-Yan 
production. 

Note that the non-uniformity of the power ex- 
pansion before integration over x does not oc- 
cur for deep inelastic scattering (DIS) processes. 
Given the correspondence of IR renormalons with 
power ultraviolet divergences of higher-twist op- 
erators, the difference between fragmentation and 
DIS must be sought in the renormalization prop- 
erties of multi-parton correlation functions that 
appear in the light-cone expansion of || . For the 
longitudinal structure function in DIS, we find 
that the quadratic power divergence at one-loop 
of the multi-parton operator 
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that is, the form of a convolution with the leading 
twist contribution. The important point to notice 
is that the operator spreads only a finite distance 
on the light-cone under renormalization. In con- 
trast, the multi-parton correlations that appear 
in fragmentation spread over the entire light-cone 
under renormalization. When the energy fraction 
x approaches zero, the operator becomes sensitive 
to very large longitudinal distances and to how 
fast the gauge fields decrease at infinity. It is this 
sensitivity to the behaviour at infinity that causes 
a 1/Q correction in the longitudinal cross section 
upon integration over x. We will return to this 
point in detail in a future publication j5j. 
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